Critical Dynamics of Thermal Conductivity at the Normal-Superconducting 

Transition 



O 

o 
o 

o 

Q 

(N 



o 
o 

Oh' 
o3 



i 

G 

O 

o 



oo 

(N 

o 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



Smitha Vishveshwara 1 and Matthew P. A. Fisher 2 
1 Department of Physics, University of California, Santa Barbara, CA 93106 
2 Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4-030 

(February 1, 2008) 

We study the effect of thermal fluctuations on the critical dynamics of the three-dimensional 
disordered superconductor across the normal-superconducting transition. We employ a phenomeno- 
logical hydrodynamic approach, and in particular find the thermal conductivity to be smooth and 
non-singular at the transition. 



I. INTRODUCTION 

Some fifteen years after the discovery of the cuprate su- 
perconductors, many tantalizing puzzles remain, particu- 
larly in regard to the ground state properties and possible 
quantum criticality as the doping is varied. It has been 
difficult to disentangle consistent analyses faithful to ex- 
periment from the zoo of microscopic theories and mod- 
els. In contrast, there have been notable successes in un- 
derstanding the phenomenology of the finite temperature 
transition into the superconducting phase, where quan- 
tum fluctuations are unimportant and progress can be 
made employing classical Landau-Ginzburg approaches. 
In particular, at optimal doping in YBCO, the observed 
critical behavior in the electrical conductivity, specific 
heat, penetration depth, and other quantities JIJ appears 
to be consistent with theoretical expectations. 

Recent attention has focused on thermal transport ex- 
periments 0, which have shed light on the low temper- 
ature transport of quasiparticles in the superconducting 
phase, and have revealed low temperature violations of 
the Wiedemann-Franz law in the normal state of the 
electron doped material suggestive of a non-Fermi liquid 
ground state ||. Here, we revisit the theory of ther- 
mal conductivity, focusing on the critical behavior near 
the finite temperature superconducting transition where 
progress is possible without the need for a microscopic 
quantum model. Older works within a BCS framework 
generalizing the Aslamazov-Larkin calculations to ther- 
mal conductivity Q| , have predicted a diverging thermal 
conductivity upon cooling into the superconductor, rem- 
iniscent of the behavior of He-4 at the A-transition. This 
appears to be at odds with experiment in the cuprates, 
which typically show a finite and rather smooth thermal 
conductivity as one cools through T c , with a large growth 
upon further cooling usually ascribed to quasi-ballistic 
quasiparticle transport. 

Our study focuses on the three-dimensional disordered 
superconductor, most appropriate to optimally doped 
YBCO which is the least two-dimensional cuprate. We 
follow the phenomenological hydrodynamic approach to 
critical dynamics pioneered by Halperin and Hohenburg 
||. Indeed, one of the early successes of this dynamical 



scaling approach was the correct description of the di- 
verging thermal conductivity near the A-transition in Hc- 
4. Here, we modify this theory to account for impurities 
and long-ranged Coulomb forces appropriate to the su- 
perconductor. Our central conclusion is that rather than 
a divergent thermal conductivity as in He-4, the thermal 
conductivity is predicted to be finite and exhibit smooth 
non-singular behavior at the superconducting transition. 
Surprisingly, the critical singularities which dominate the 
electrical conductivity and other thermodynamic proper- 
ties are found to be completely absent from the thermal 
conductivity. There is effectively a decoupling between 
the thermal and electrical transport coefficients upon ap- 
proaching the superconducting phase, in strong contrast 
to the universal Lorenz ratio relating these two transport 
coefficients in a conventional metal. 



II. THERMAL CONDUCTIVITY 

The thermal conductivity k, relates the heat current Q 
to an applied thermal gradient VT under the condition 
of no particle flow, 



Q = -kVT, j = 0, 



(1) 



where j is the particle current. Within linear response 
one has an Einstein relation for the thermal conductivity, 



DtCv, 



(2) 



where Dt is the thermal diffusion constant, and Cy the 
specific heat at constant volume. 

In systems with a condensed ground state such as 
superfluid He-4, the situation is a little more involved 
since energy can couple to the order parameter in a non- 
dissipative fashion. This leads to a ballistic wave-like pro- 
pogation of heat - second sound - in the superfluid state 
, and an infinite thermal conductivity. In the normal 
state, while heat does propogate diffusively, its associ- 
ated thermal conductivity diverges due to critical fluctu- 
ations upon approaching the transition temperature. As 
we shall discuss, in contrast to its sister system of super- 
fluid He-4, the impure superconductor can only support 
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a diffusive rather than an oscillatory heat mode, and the 
thermal conductivity k is finite at all temperatures. 

In the high T c superconductors, the role of thermal 
fluctuations near the superconducting transition is dras- 
tically enhanced relative to their low T c counterparts. 
Being generally anisotropic, and strongly type II, the 
Ginzburg criterion shows that critical fluctuations are 
present over a relatively wide range of temperatures, per- 
haps as large as 5-10 Kelvin M. Within this temperature 
window we can dispense with microscopic models, and 
appeal to a critical hydrodynamic approach. 



III. SUPERFLUID HYDRODYNAMICS 
REVISITED 

To identify the appropriate model near criticality, we 
follow Ref. [^), and first focus on the hydrodynamics in 
the normal state. Since the frequency u> of the 'slow' hy- 
drodynamic modes vanishes with wave- vector q, one need 
only study the dynamics of conserved densities which 
cannot relax rapidly on long lengthscales. We consider 
first an uncharged and pure fluid, such as He-4. In this 
case there are five such conserved densities: three com- 
ponents of the particle current density (J) and the energy 
(e), in addition to the particle number density (p). Asso- 
ciated with these five conserved densities are five hydro- 
dynamic modes: propagating first sound which involves 
the density and longitudinal current density (counting 
as two modes with lu = icq) and three diffusive modes 
(with u> = iDq 2 ) - energy and two transverse current den- 
sity modes. To see this, consider the continuity equation, 



d t e 



D T V 2 e 



(6) 



dp ^ 



0, 



(3) 



and the Navier Stokes equation for current density con- 
servation linearized for small current density: 



% + -Vp = vV 2 ], 
at m 



(4) 



where p is the pressure , m the mass of the constituent 
He-4 particle, and v the kinematic viscosity. One can also 
define the fluid velocity field v, which satisfies j = pv. 
Upon taking the divergence of EqnJJ, combining with 
the continuity equation and linearizing for small current 
density and density variations, one obtains first sound 
l^js. Eqn.[| also implies two dif- 



with velocity, c = 
fusive transverse current density modes. In addition, the 
energy density e is conserved: 



d( 

^ + V.Q = 0, 



(5) 



where Q is the heat current. Upon combining Eqn.|^ with 
the expression that defines the thermal conductivity, one 
obtains the heat diffusion equation, 



where we have used Cy = de/dT and the Einstein rela- 
tion relating k and the thermal diffusion coefficient, Dt- 
On approaching the A-transition into the supcrfluid, 
one must augment this hydrodynamics with the slow dy- 
namical relaxation of the order parameter. This leads 
to six hydrodynamic modes in the superfluid phase: first 
sound, second sound (involving the order parameter and 
predominantly energy) and the two diffusive transverse 
current density modes. In the superfluid one adopts a 
"two- fluid" description, in which the total density is de- 
composed into a superfluid and a normal fluid compo- 
nent: p = p s + p n . Moreover, one introduces a second 
velocity field - the superfluid velocity v s , in addition to 
the "normal" fluid velocity denoted v n , and the total 
current density is decomposed as j — p s v s + p n v n . The 
superfluid velocity field is taken to satisfy a Josephson 
relation, 



dv s -» 
m— + V/i = 0, 
at 



(7) 



where p is the chemical potential [|J . Diffusion of energy 
in the normal fluid is replaced by second sound oscilla- 
tions of normal and superfluid. As the superfluid carries 
no entropy, the oscillations involve an entropy wave asso- 
ciated with the normal fluid. This is easiest seen in the 
absence of dissipation {y = 0, K = 0), where conservation 
of entropy density s, maybe expressed in terms of the as- 
sociated entropy current sv n carried by the normal fluid 
as 



ds _ 



0. 



(8) 



One can combine Eq.||,[I[f7]||, and use the fact that local 
changes in chemical potential are related to local thermal 
gradients via the thermodynamic relation 



pdp = —sdT + dp. 



(9) 



One then finds the superfluid velocity and the entropy 
participating in second sound motion with velocity c s = 



j'... .s- 



First sound and the two diffusive transverse current 
density modes are also present in the superfluid phase, 
giving the anticipated six hydrodynamic modes. 

Turning next to the critical dynamics at the A- 
transition, one needs to take into account purely those 
hydrodynamic modes which go soft. Coming from the 
superfluid side, the only such mode is the second-sound 
mode since its velocity vanishes along with the superfluid 
density. The first sound velocity remains finite through 
the transition, so that first sound oscillates much more 
quickly than second sound at the same wavevector. One 
thereby argues that first sound should not enter into the 
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low frequency critical dynamics. Moreover, couplings of 
transverse current density modes are found to be irrel- 
evant Thus for the critical dynamics, one can fo- 
cus on the conserved heat density and the order param- 
eter, which together comprise second sound in the su- 
pcrfluid. More precisely, the non-dissipative coupling of 
energy and the order parameter below T c , which results 
in second sound, is captured by the Poisson relationship 



(10) 



where M is the net combination of the energy and par- 
ticle number densities taking part in the oscillation, and 
tjj is the complex order parameter. The critical dynam- 
ics of such a system is described by a generalized time- 
dependent Ginzburg-Landau type model involving these 
two fields, and was denoted as Model F in Ref. || . 



IV. COULOMB INTERACTIONS AND 
IMPURITIES 

Turning to the case of the charged, pure fluid, which 
is appropriate in the normal phase for a pure metal, the 
presence of an electric field E gives rise to acceleration. 
Thus, the r.h.s of Eq.|| acquires a term j^pE, and that 
of Eqjl] acquires a term eE, where 'e' is the charge of 
the electron. In addition, distortions of the charged fluid 
itself give rise to an electric field described by Poisson's 
equation, 



\7.E = AneSp, 



(11) 



where 5p corresponds to density distortions. Conse- 
quently, in three dimensions, the sound mode involving 
density and longitudinal current density becomes a high- 



energy plasmon with frequency uj p = y 47r ^ 2p in the limit 
q — > 0, and thus drops out of the hydrodynamic descrip- 
tion. The transverse diffusive current density modes re- 
main. Importantly, even with the modifications associ- 
ated with the presence of charge, one finds the diffusive 
heat mode above T c . 

Below T c , in the pure three dimensional superconduc- 
tor, the low frequency, long wavelength second sound 
mode survives because it involves no net density distor- 
tions. To see this, we note that due to the gapped plas- 
mon the density is not a hydrodynamic variable, so one 
can effectively set bp — 0. Upon linearizing the continu- 
ity equation EqJ^ for small velocities, this implies that 
V • v n = — (p s //0„)V • v s . Upon inserting this into the lin- 
earized entropy conservation Eqn.^J, and combining with 
the divergence of the Josephson relation Eq.0, one finds 
that 



— + PsS y 2 
dt 2 p n m 



Finally, using the thermodynamic relation Eqn.[| one ar- 
rives at the wave equation for the entropy density pro- 
pogating once more with velocity c s — * / p °fi, , . In 



P m Pn(f§0 ' 

addition to second sound, one expects the two diffusive 
transverse current density modes to be present in the 
charged superfluid, just as in He-4. 

Thus, essentially the sole effect of Coulomb interac- 
tions on the superfluid hydrodynamics is the conversion 
of first sound into the non-hydrodynamic plasmon mode. 
Since first sound was argued in any event to decouple 
from the critical dynamics at the A-transition, the criti- 
cal dynamics in the charged superfluid is expected to be 
described by the same theory - that is model F, except 
with the density M referring to pure energy. 

Consider next the case of an uncharged fluid in the 
presence of impurities (eg. He-4 absorbed in a porous 
medium). Impurities violate momentum conservation 
and this leads to a dramatic modification of the hydrody- 
namics. In particular, with only density and energy being 
conserved, one expects two hydrodynamic modes above 
T c and three below. Absence of momentum conservation 
can be explicitly captured by 



dj n 
dt 



Jn 
T 



(13) 



0. 



(12) 



where 1/r is the decay rate of current density j n . Then, 
above T c , the first sound mode involving longitudinal cur- 
rent density and density cannot propagate at low frequen- 
cies. It is replaced by a damped current density mode 
(lu = —i/r) and a diffusive density mode with diffusion 
constant D p — c 2 t, where 'c' is the velocity of first sound 
in the pure case. Energy continues to be a diffusive mode 
in the presence of impurities. 

Below T c , the uncharged impure system once again has 
a superfluid component moving with velocity v s . The su- 
perfluid can couple to density in a non-dissipative fash- 
ion, but not to entropy, which is not even conserved, and 
is still mainly associated with the normal fluid that can 
no longer propagate ballistically. We can thus simply 
drop the normal fluid velocity from the continuity equa- 
tion EqnJ^, which upon linearization and combination 
with the Josephson relation gives, dfp = (p s /m)V 2 /i. 
This describes a fourth sound mode involving purely su- 
perfluid oscillations propagating with velocity Ci rnp = 
I p g p . In addition, the diffusive energy mode persists 

below T c . 

To model the critical dynamics one needs to re- 
interpret the conserved density M in model F as the 
particle number density, p, and then augment the model 
with an additional conserved energy density which is dif- 
fusive both above and below T c . The conserved energy 
density will be coupled to the order parameter -0 via a 
term in the free energy of the form e|f/ , | 2 - This will lead 
to a coupling in the resulting time-dependent equations 
of motion. 
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Finally, we arrive at the case of interest -the impure 
superconductor-, which we model by the charged, impure 
superfluid. The hydrodynamic modes may be obtained 
by incorporating the modifications described above for 
charge and impurities to the linearized hydrodynamic 
equations of the He-4 system. The normal state describes 
an impure metal. Due to the long-ranged Coulomb inter- 
actions the density, although conserved, is gapped up at 
the plasma frequency and thus again drops out of the hy- 
drodynamic description. Moreover, with impurities, the 
transverse current density modes are also damped. So 
above T c , in this system of the impure metal, associated 
with the only remaining hydrodynamic variable one has 
a diffusive thermal mode. 

Below T c , superconducting order gives rise to a second 
mode. As in the case of the uncharged, impure fluid, 
the order parameter cannot couple non-dissipatively to 
heat. But unlike in this case, it cannot even do so with 
density, which as in the pure charged case, is no longer 
a hydrodynamic variable since its distortions cost elec- 
trostatic energy. As a result, there exist two diffusive 
modes related to the energy and order parameter, and 
no oscillatory hydrodynamic modes 1 . 

While our hydrodynamic treatment employed simple 
linearized equations appropriate for mode analysis, an 
exhaustive study along the lines of the case of He-4 
may prove useful for the case of the impure superconduc- 
tor. 



V. CRITIC ALITY 

We are now equipped to model the critical dynamics 
of the dirty superconductor. We have seen that above 
T c , there only exists a diffusive energy mode. Below 
T c , in striking contrast to superfluid He A , instead of a 
second sound oscillatory mode, there exists one diffu- 
sive mode associated with conserved energy density, and 
one diffusive mode associated with the order parameter, 
which takes the form of a complex scalar. The simplest 
phenomcnological model incorporating these ingredients, 
and all possible relevant couplings, is described by Model 
C of Ref. H , and is analyzed in detail for its critical dy- 
namics in Ref. [jn]]. ft is defined by the following set of 
equations of motion involving the complex order param- 
eter ip and the conserved energy density: 



dt 



SF 



(14) 



de 
dt 



5e 



(15) 



F„ = / d d x{K \^\ 2 + Uo H 4 + i|Wf + 7 oM 2 f 



1 



1,2s 



(16) 



Here Tq and kq are the bare transport co-efficients for the 
order parameter and the energy, respectively. An exter- 
nal source field, 5 has been included in Eqn.[l6|, and rj and 
£ are Langevin noise sources appropriate for ip (complex) 
and e (real), respectively. The energy density and the su- 
perconducting order parameter are coupled together via 
the coupling constant 70. 

In equilibrium, i\) and e minimize the functional Fq. In 
a functional integral formulation applicable for the equi- 
librium distribution fll|| , one can integrate over the en- 
ergy field, e. The resulting functional of the order param- 
eter i\) displays the statics of 3D XY critical behaviour, 
appropriate for an extreme type II superconductor (see 
for e.g. Ref ||J^]). The Harris criterion (see for e.g. Ref 
|l2|) shows that this holds true even in the presence of 
disorder. Thus, the pure 3D XY model provides an ap- 
propriate description for critical static properties of the 
dirty superconductor. 

We next turn to the critical dynamics. To obtain the 
thermal conductivity, one first defines the energy density 
linear response function x(cf, lo), by the relation 



< e(<7, > s = x(<f, u))S(q, uj). 



(17) 



Quite generally, in such a diffusive system, one expects 
this response function at low frequencies and wavevectors 
to take the standard form: 



X(<2>) = 



— iuj + Dxq 2 ' 



(18) 



where k is the thermal conductivity and Dt the thermal 
diffusion constant. Thus, the thermal conductivity can 
be extracted as, 



lim q 

q-^0 









u=0- 



(19) 



In the absence of coupling to the order parameter, \ 
can be readily computed from the (linear) equation of 
motion for e, which gives, 



-iu> 



noq 



(20) 



1 Oscillatory modes involving normal and superfluid densities 
were observed in aluminium films ||. However, modified hy- 
drodynamic treatments Q show that except for a narrow win- 
dow of momenta, these Carlson-Goldman modes are damped. 
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which has the general diffusive form as in Eqn.[18] with 
the identifications k = kq and Dt — kq/Cq. 

The effect of the fluctuating order parameter on the 
thermal conductivity k, can be studied by treating the 
couplings uq and 70 perturbatively. Specifically, one can 
set up a perturbative expansion for the "self-energy" (£) 
which modifies the thermal response function: 



NSF Grants DMR-97-04005, DMR95-28578 and PHY94- 
07194. 



X 



(21) 



As is clear from Eqn.[l9[ a renormalization of the bare 
thermal conductivity kq requires a contribution to the 
self energy of the form X ~ —iuj/q 2 , which is divergent 
as q — > for fixed non-zero to. But as argued by Ref. 
flEl), away from criticality, with tq > 0, the self energy 
will be finite in the q — > limit at all orders in perturba- 
tion theory, being protected in the infrared by r and in 
the ultraviolet by a high momentum cutoff. Thus, these 
perturbations can generate no corrections to k, and the 
thermal conductivity will be given exactly by the " bare" 
parameter n a 11 1. Since n is a coupling constant which 



depends on short distance physics only, it will necessar- 
ily be a smooth function of temperature. This thereby 
establishes that the thermal conductivity should be non- 
singular and smooth as one cools through the supercon- 
ducting transition. 

It is worth emphasizing that although the thermal con- 
ductivity itself is non-singular at the transition, critical 
singularities driven by the order parameter flucutations 
will generally enter into the dynamical relaxation of the 
energy. Smoothness of the thermal conductivity is inti- 
mately tied to the fact that it is a zero-frequency, zero- 
wavevector quantity. Indeed, the dependence of n on 
wavevector or frequency is expected to be singular at the 
critical point. 

To recapitulate and emphasize, fluctuations in the su- 
perconducting order parameter will certainly affect the 
time dependent relaxation of the energy. However, the 
thermal conductivity - the transport co-efficient associ- 
ated with energy propagation - remains smooth and finite 
across the transition. 
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